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We study the extended Hubbard model with both on-site {U) and nearest neighbor {V) Coulomb 
repulsion using the exact diagonalization method within the dynamical mean field theory. For a 
fixed U {U — 2.0), the T — n phase-diagrams are obtained for V = 1.4 and V = 1.2, at which 
the ground state of n = 1/2 system is charge-ordered and charge-disordered, respectively. In both 
cases, robust charge order is found at finite temperature and in an extended filling regime around 
n = 1/2. The order parameter changes non-monotonously with temperature. For V = 1.4, phase 
separation between charge-ordered and charge-disordered phases is observed in the low temperature 
and n < 0.5 regime. It is described by an "S"-shaped structure of the n — /i curve. For V — 1.2, 
the ground state is charge-disordered, and a reentrant charge-ordering transition is observed for 
0.42 < n < 0.68. Relevance of our results to experiments for doped manganites is discussed. 



I. INTRODUCTION 

Charge-ordering is a fascinating research topic in con- 
densed matter physics. In recent years, charge or- 
der and the related spin and orbital order in doped 
manganites have attracted much interest. In the sys- 
tem Lai-xCuxMnOs {x > 0.5), charge order in the 
ground state is enhanced as x increases. The correspond- 
ing transition temperature Tco is higher than the Neel 
temperature^, and there are strong charge-order fluctu- 
ations at temperatures above Tco (Rcf.^). In particular, 
a reentrant charge-ordering transition has been observed 
in systems such as /'ro.65(Cao.7»S'''o.3)o.35-^'^03 (Ref.'^) 
and LaSr^Mn^Oj (Ref.*^^). The origin of these charge- 
ordering transitions lies in the complex interplay between 
orbital and lattice degrees of freedom in the manganites, 
and is presently under intensive study. In other systems 
such as the heavy fermion system Yh^Asj, (Ref.^), the 
quasi-one dimensional material NaV20c, (Ref.^'*), and 
the superconducting layered organic molecular crystal 
K-{BEDT-TTF)2X (Ref.^), charge-ordering is closely 
related to the specific properties of the system. 

The charge ordering in the above stated systems has 
different physical origins, and cannot be explained within 
a single theory. From a theoretical point of view, an obvi- 
ous cause of charge-ordering is the short-range Coulomb 
repulsion between electrons. The simplest model that in- 
cludes this interaction is the extended Hubbard model 
which contains the kinetic term and the on-site and 
nearest-neighbor Coulomb repulsion. Despite the sim- 
plicity of this model, recent studies found that it can 
explain some characteristics of the experimental obser- 
vations in doped manganites^*^"^^. These studies re- 
vealed that many interesting effects arise simply from 



pure Coulomb repulsion and charge fluctuations. In 
particular, using dynamical mean-field theory (DMFT), 
Pietig et al. found that the quarter-filled extended Hub- 
bard model exhibits a reentrant charge-ordering tran- 
sition near a critical value Vc of the nearest-neighbor 
repulsion^ ^ (Fig. 1). It coincides with what was observed 
in the doped manganites Pro. 65(C'ao.7<S'ro. 3)0.35^/7103 
(Ref.3) and LaSr2Mn207 (Ref.^'^). This indicates that 
some properties of charge order may be independent of 
the concrete microscopic mechanism, and hence can be 
studied using simplified models such as the quarter-filled 
extended Hubbard model. Along this line, other authors 
also studied this problem^'^^^^ . In lower dimensions, the 
model has also been extensively studied in various con- 
texts using different methods^^. 

Experimentally, in the doped manganites, charge- 
ordering was observed not only at x = 0.5, but also 
in a broad doping regime. Also the reentrant behav- 
ior was observed in systems away from quarter-filling. 
In fact, Pr(iQ5{CaojSro,3)o,35Mn03 has an electron fill- 
ing oi n = 0.65 > 0.5. Another feature of the experi- 
mental observation is the inhomogeneous coexistence of 
charge-ordered and charge-disordered phases. In the sys- 
tem Ri-xCaxMnOz (R = La, Nd, Bi etc.), many ex- 
periments showed that phase separation (PS) between a 
ferromagnetic charge-disordered phase and an antiferro- 
magnetic charge-ordered phase exists at dopings rang- 
ing from X = 0.33 (Ref.^*) to x = 0.82 (Ref.^^). Al- 
though detailed explanations should take into account 
the spin and orbital degrees of freedom, simplified mod- 
els also give remarkably similar results. In the spin- 
less fermion model, when electrons (holes) are added to 
a half-filled system (n = 1/2), an electron (hole) rich 
charge-disordered phase is separated from the charge- 
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ordered background^". 

In this paper, motivated by the experimental observa- 
tions as well as by the work of Pietig et al.^^, we study 
the properties of a single-band extended Hubbard model 
away from but near quarter-filling. We confine our study 
to the two-sublattice case and do not consider a possi- 
ble incommensurate ordering. We are interested in the 
following two issues of charge-ordering. First, the stabil- 
ity of the reentrant charge-order transition in the whole 
doping regime. Second, the possibility of PS between 
charge-disordered and charge-ordered phases. These is- 
sues are relevant to the charge-ordering and PS phenom- 
ena observed in doped manganites. Theoretical studies 
in the past focused only on fillings at or near n = 1.0 and 
n = 0.5. Despite the intensive studies, these two issues 
remain unclear. 
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FIG. 1. Schematic T — V phase diagram for the extended 
Hubbard model at C/ = 2.0 and n = 1/2 (taken from Ref.^^, 
with V rcscalcd.). CO and CD denote charge- ordered and 
charge- disordered, respectively. The vertical dashed lines 
mark out the positions of T — n planes studied in this pam- 
per, i.e., V = 1.2 and V = 1.4, respectively. 

Using DMFT together with the exact diagonalization 
technique, we study the paramagnetic phase diagram in 
two T — n planes for a fixed [/-value U = 2. As shown 
in Fig. 1, these two planes cross the T — V plane at 
V = 1.2 and V = 1.4, which are larger and smaller than 
Vc, respectively (Vc is the critical value of the nearest- 
neighbor repulsion for a quarter-filled system, which 
separates charge-ordered and charge-disordered ground 
states. Vo « 1.32W^ (Ref"), see Fig. 1.). We found 
that a reentrant charge-ordering transition exists in an 
extended regime of the electron density n near quarter- 
filling. 

For V > Vc, the order-paramter of charge-ordering 
changes non-monotonously as temperature decreases. 
The ground state is still charge-ordered for n > 0.5 
regime and charge-disordered for n < 0.5. In the for- 
mer regime. We also find the PS between charge-ordered 
and charge-disordered phases. 



For V a little smaller than Vc, the ground state is 
charge-disordered for any n. The charge ordering ex- 
ists only at finite temperatures and is most robust in the 
n > 0.5 regime. At the lower critical temperature of 
the reentrant transition, the order-parameter disappears 
rather abruptly. Phase diagrams in these two planes are 
plotted. 

In Sec. II, we describe the model and the method used 
in this work. In Sec. Ill, our results for V = lAW > Vc 
are presented. Phase separation between charge-ordered 
and charge- disordered phases are discussed. Their ther- 
modynamical structure is compared with other kinds of 
PS's. Our results for V = 1.2W < Vc are presented in 
Sec. IV. A summary is given in Sec. V. 

II. MODEL AND METHOD 

The Hamiltonian of the single-band extended Hubbard 
model has the form 

+ ^rii^j - yu^nj. (1) 

In Eq.(l), J2{i j) indicates the sum over nearest-neighbor 
sites i and j independently. J2{ij) indicates the sum 
over nearest pairs. Hence there is the relation = 
^/2J2{ij)- U and V are on-site and nearest-neighbor 
Coulomb repulsion, respectively, and ji is the chemical 
potential. Here we use a Bethe lattice which produces 
a semicircular density of states for free electrons in the 
limit of large coordination number. In this limit, the 
inter-site coupling terms in Eq.(l) are replaced by the 
corresponding Hartree term. After a proper rescaling of 
V: V ^V/Z (Here Z is the number of nearest neighbors) 
and neglecting the constant term, we obtain the following 
mean-field Hamiltonian: 

-^(M-l^(ni+5))ni, (2) 

i 

where i -I- (5 is the nearest site of i. Note that our scaling 
is different from that of Ref.^^, where the scaling V 
2V/Z is used. So in this paper the value of V is twice as 
large as the corresponding one in Ref.^^. 

To describe the charge-ordered phase, we divide the 
Bethe lattice into two sublattices. Correspondingly, 
within DMFT, the model Eq.(2) is mapped onto two un- 
coupled effective Anderson impurity models, 

Ns-l 

-f^A.imp = ^ \ihka\^^aAka +VKk{a\ka'^ha + h.C.) 
k=l,a 

+ Un'^nl-{n-V{n^))nl. (3) 
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Here, A = A,B refers to the two sublattices {A = B, 
B = A) and {eAfe, VXfc} are effective parameters describ- 
ing the bath. They are related to the Weiss function 
^AO (*^") through 
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map 



Ak 



CAfe 



(4) 



We use the full exact diagonalization method to calcu- 
late the impurity Green's function for this model. The 
number of sites Ng = 5 and Ns =6 are found to be suffi- 
cient for the calculations in this paper. The free density 
of states is given by: 



D{e) 



VW^ - e2 {\e\<W). 



(5) 



We set = 1 as the energy unit. The DMFT self- 
consistency equations for the Bethe lattice are the given 
by 

Sao {i^n)dys =iuJn + IJ'-V (n-x) —Gj;{iu>n), (6) 

where A = A, B. Equations (3), (4), and (6) form a 
set of closed self-consistent DMFT equations. After the 
Green's function and Weiss function arc calculated, the 
new set of effective parameters for the impurity model is 
obtained via a minimization procedure^^ : 
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"max 

E 

n=0,A 



OAo{i^n)map — QAo{i^n)dys\ ■ (7) 



The self-consistent equations are solved iteratively. Af- 
ter the iteration converges, we calculate the electron den- 
sities on the two sublattices, ua and n^, which in turn 
give the average electron density n = {ua + jt-b) /2 and 
the order parameter \nA — hbI- The contributions to the 
total energy per lattice site E — Et + Eu + Ey from dif- 
ferent parts of the Hamiltonian are also calculated. The 
kinetic energy Et, on-site repulsion energy Eu, and the 
inter-site repulsion energy Ey are given by: 



(8) 



n,A 



where 

(«w„) = ia;„ -I- - g~^Q (ia;„) (A = A,B), (9) 
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Ev = -^nATlB ■ 



(11) 



Phase separation has been studied extensively for mod- 
els of strongly correlated electrons, such as the Hubbard 
model22-24^ the t - J modeP^'^^, the Falicov-Kimball 
modcl^^and the double-exchange modcP*^"''". A stan- 
dard criterion for the PS is the discontinuous jump in 
the curve of electron density n versus chemical poten- 
tial yLt (Ref.^^). However, in the DMFT study of double- 
exchange systems, we have observed PS also through the 
multiple- valued structure in a continuous n — ^ curve^^. 
As described in the next section, we find in a certain 
temperature regime for F = 1.4 an "S" or "Z" -shaped 
multiple- valued structure in the curves of ua, tlb and n 
versus /U. Such curves contain the full information about 
the PS, including the metastable phase and the first-order 
phase transition line. To do this, we introduce^^ a self- 
consistently determined quantity /i', 



fjL = jjL — \ {nA — A) . 



(12) 



where A and A are tunable parameters. The value of 
Ha is dependent on the chemical potential fi through 
Ua = F{ii) and the functional dependence is de- 

termined by the DMFT calculations. For a given /z and 
in each DMFT iteration, we first calculate /i' and use /i' 
instead of ji in the ordinary DMFT scheme to produce 
the local Green's function, and then extract the CAk and 
VXfe- After that, the new value of ua is calculated. The 
iteration is carried on until convergence is reached. This 
is equivalent to simultaneously solve the DMFT equation 
and the following equation. 



nA = F[^l-\{nA- A)]. 



(13) 



If the function F{ijl) has a multiple-valued regime, by 
selecting appropriate parameters A and A, the ua — IJ- 
curve self-consistently determined by Eq. (13) may be- 
come single-valued. In this way, our calculation avoids 
the numerical instabilities induced by the multiple- valued 
structure of n — /x curve. For each ^, we first solve Eq. 
(13) together with the DMFT equations, then calculate 
the thermodynamic quantities Q that we are interested 
in. After the data is obtained for each fi, we plot the 
quantities Q with respect to the argument /z' to recover 
the physical curves Q ~ Fq (p') that correspond to the 
orginal Hamiltonian. The final results should be inde- 
pendent of the parameters A and A, if only they are in 
an appropriate regime. In the framework of DMFT, this 
transformation scheme has been used in the study of PS 
in the double exchange modeP^ and of the Mott-Hubbard 
transition in the Hubbard modeP^. 
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[S-aJ (iw„) - G^^ {iujn)] Ga (iw„) e 
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(10) 



III. RESULTS AND DISCUSSION 

A. Phase separation: V = 1.4 > Vc 

In this section, we discuss our results for V = 1.4, 
which is a little larger than the zero-temperature crit- 
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ical value Vc « 1.32 (See Fig.l). For this interaction 
strength, as is shown in Fig. 1, the exactly quarter-filled 
system has a charge-ordered ground state. The charge 
order persists up to T « 0.26, and no reentrant tran- 
sition was found. In order to study the system away 
from quarter filling, at a lower temperature T = 0.05, 
we change the chemical potential and calculate charge 
densities. Fig. 2(a) shows our result for the sublattice 
charge densities ua, ub, and average charge density n 
as functions of /x. As we have expected, charge order 
exists in some finite regime of doping around n = 1/2: 
0.44 < n < 0.73. This regime is not symmetric about 
n = 1/2, since the model (1) does not have particle-hole 
symmetry at this point. This differs from the spinless 
fermion model^". Here, the charge-ordered regime ex- 
tends to larger densities, where the inter-site Coulomb re- 
pulsion is more effective. However, the value of — nsl 
is largest at n = 1/2. 
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FIG. 2. (a) electron density n versus chemical potential fj, 
for U=2.0, V=1.4 and T=0.05. The squares, dots, and dia- 
monds are for ua-, ub, and n, respectively. Calculation results 
with Ns=6 (hollow symbols) and Ns=5 (cross-filled symbols) 
agree well. Linos arc for guiding eyes, (b) Contributions to 
the energy per lattice site as functions of n: Et, Ey, Eu and 
E denote kinetic energy, nerest neighbor repulsion energy, 
on-site repulsion energy, and the total energy, respectively. 
The regime between the dotted vertical lines is unstable to- 
wards phase separation. Charge order exists between the solid 
vertical lines. 



One dominant feature of Fig. 2(a) is the "S"-shaped 
multiple- valued structure in the n — fi, ua — /i and H 
curves. This is direct evidence for PS in the extended 



Hubbard model near quarter filling. When n increases, 
the system first goes from a charge-disordered phase into 
a charge-ordered phase through a first-order transition, 
which occurs in the regime n = 0.4~0.5. Asn increases 
further, the charge order disappears continuously at n « 
0.73. The results obtained using Ng = 5 and Ng = 6 
agree very well in most part of the curve. There is only a 
slight deviation near the second-order transition point at 
n = 0.44. This indicates that for this problem, iV^ = 5 is 
sufficient to obtain qualitatively correct conclusions. In 
the following, all our results are from = 5 calculations. 

In the multiple- valued regime of the n — n curve, the 
DMFT self-consistency equations have three solutions for 
a fixed fi. We use the thcrmodynamical grand potential 
LO (T, jj) to evaluate the relative stability of these solu- 
tions: 



(T, m) = w (T, ^o) 



n (T, /i') dfi' . 



(14) 



The actual first-order transition point fic is determined 
by a Maxwell construction, i.e., by solving the equation 
Wi {T,fic) = u)2 {T,Hc)- At fi = /Xc, two phases coexist, 
with their respective volumes determined by the nominal 
electron density of the system. At T = 0.05, one of the 
coexisting phases is charge-disordered with ni = 0.415, 
the other one is charge-ordered with n2 = 0.486. The 
third solution with intermediate n is charge-ordered, but 
due to its negative compressibility and highest grand po- 
tential, it is unstable with respect to the others. For /i 
away from but near fi^, there is only one stable phase in 
the system, either ordered or disordered. The other two 
solutions have higher grand potential and are metastable. 
These metastable phases may be detected by hysteresis 
experiments. 

n = 0.44 is the lower critical density at which charge 
ordering occurs. At this point, the n ^ ii curve turns 
backwards sharply, and n = — 9a;(T, /z) /djj is continu- 
ous while dn/dfj, is discontinuous. Therefore it is iden- 
tified as a second-order transition point between charge- 
disordered and charge-ordered phases. In this respect, 
it is same as the n = 0.73 point. However, due to the 
higher grand potential, the second-order charge ordering 
transition only exists at metastable level. To understand 
the PS better, we show the n-dependence of the energy 
in Fig. 2(b). The on-site repulsion energy Eu has a small 
contribution. It does not change much at the boundary 
to the charge ordering regime. This is an indication that 
for U = 2, double occupancy is small in the doping range 
around n ~ 0.5, and that it plays a minor role for PS. In 
contrast, the kinetic energy Et and inter-site repulsion 
energy Ey are more sensitive to the charge density and 
long-range charge ordering. If there is no charge ordering, 
Ey behaves as Ey oc n^, while Et recovers the behavior 
of the Hubbard model; this means that with increasing 
n, it decreases in the low filling area and increases near 
half filling. 

In Fig. 2(b), it is seen that the charge order in the 
regime 0.44 < n < 0.73 strongly reduces Ey while it 
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increases Et with respect to their value in disordered 
phase. The charge ordering transition at n = 0.44 causes 
drastic changes in Ey and Et, which arc then natu- 
rally related to the occurence of PS. In contrast, at the 
other transition point n = 0.73, the energies change 
more smoothly. Because the opposite contributions from 
Ev and Et almost cancel, the total energy E increases 
monotonously and smoothly with increasing n. However, 
the quick change of Ey and Et near n = 0.44 and their 
competition leads to a small convex part in the E — n 
curve, where d'^Ejdr? < 0. When taking into account 
the entropy contribution to the free energy, the convex 
structure will be enforced by a small amount TS ~ 10~^. 
Therefore, at the temperature T = 0.05, the E — T curve 
already represents the behavior of free energy F (T, n) . 
The convex structure observed here is therefore consis- 
tent with PS obtained using the n — /i criterion. From 
Fig. 2(b), it is clear that PS is closely related to the 
charge ordering transition. It directly results from the 
charge-ordering-induced strong competition between Ey 
and Et- 
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FIG. 3. Electron density n versus chemical potential for 
U = 2.0, V = 1.4, and (a) T = 0.01, (b) T = 0.1, (c) T = 0.13. 
The upper and lower dotted lines and the solid line are for 
TiA, riB, and n, respectively. 

In the following, we study the temperature dependence 
of charge order and PS. In Fig. 3, three n — ^ (as well as 
UA — H and hb — n) isotherms are shown for T = 0.01, 
0.1, and 0.13. Up to temperature T = 0.13, the charge 
order is quite robust in the intermediate density regime 
from n 0.3 to n ~ 0.7. In contrast, PS is stable 



only at much lower temperatures. Compared with the 
T = 0.05 curve in Fig. 2(a), the multiple-valued struc- 
ture is more pronounced at T = 0.01, while it disappears 
at higher temperature. As temperature decreases, the 
multiple-valued structure in n — /x curve is compressed 
along the n axis, but that in the ua — M and tlb — H 
curves docs not change much. As a result, the aver- 
age electron densities of the coexisting two phases, which 
are determined through Maxwell construction, get closer 
and both move towards n = 0.5. On the other hand, 
as temperature increases, the multiple- valued part of the 
n — fj, curve shrinks along the fj, axis until it disappears at 
T w 0.1(Fig.3(b)). At this temperature, PS disappears 
and the slope of the n — /i curve diverges at the second- 
order charge ordering transition, leading to strong fluctu- 
ations in charge density as well as in the order parameter 
\nA ~ nB\- The critical end point of the PS is estimated 
to be n « 0.39, T « 0.1. 

It is interesting to note that the low temperature be- 
havior of PS is not unique to this model. Similar behavior 
has been observed in the DMFT study of first-order phase 
transitions in other strongly correlated electron mod- 
els. The ferromagnetic-paramagnetic (FM-PM) PS in the 
double-exchange model^^ and the Mott-Hubbard metal- 
insulator transition in the half-filled Hubbard modeP^ 
are both typical first-order phase transitions. Within 
DMFT, they are described by continuous "S" or "Z"- 
shaped curves that are very similar with the n — fj, 
curves presented here. As the temperature decreases, 
all these structures are compressed in the vertical direc- 
tion. The " order parameter" of the first-order transitions 
(|n_FM — npM\ for FM-PM PS, and double occupation 
difference \DMet — Dins\ for the Mott-Hubbard transi- 
tion) obtained from a Maxwell construction reduces to 
zero in the limit of T — > 0. 

Here, though the lowest temperature that we study 
is T = 0.01, the main tendency is clearly that as tem- 
perature decreases, the upper two branches of the "S"- 
shapcd curve (as shown in Fig. 3(a)) tend to merge. 
Hence we expect that as T ^ 0, the density difference 
between the two coexisting phases reduces to zero, sim- 
ilar with our previous findings in other systems. The 
two coexisting phases have the same density n = 0.5 and 
the same energy. In such a scenario, the ground state 
of the V = lA system is singular at n = 0.5. At this 
point, a charge- disordered phase and a charge-ordered 
phase, both with avarage density n = 0.5, coexist in 
random volume proportion. An infinitesimal amount of 
additional holes in the system will destroy the phase- 
separated ground state and turn it into charge-disordered 
state, while electrons will turn it into charge-ordered 
state with \nA — nsj ~ 0.7. Again, this is different from 
the spinless fermion system, for which both additional 
holes and electrons doped into the charge-ordered ground 
state at n = 0.5 cause PS^°. 

In Fig. 4, the T — n phase diagram for U = 2.0 and 
V = 1.4 is shown. It is seen that charge order is rather 
robust in the filling regime around n = 0.5. At zero 
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temperature, the charge-ordered ground state extends 
from n = 0.5 to n « 0.69. This range first expands, 
and then shrinks as temperature increases, leading to a 
reentrant transition of charge ordering in certain dop- 
ing regimes. The highest charge ordering transition tem- 
perature Tco.max ~ 0.28 is reached at n = 0.5. The 
PS area hes near the lower transition temperature line 
on the left side. When both the nominal electron filling 
and the temperature lie in this charge-disordered 
phase with lower electron density will be separated from 
a charge-ordered phase with higher electron density. The 
two boundaries of this coexisting 3XG3j^ clS shown by the 
solid lines in the figure, meet at two end points. The finite 
temperature end point is located at (Tc ~ 0.1, nc ~ 0.39), 
which is easily seen for ordinary first-order phase transi- 
tions such as liquid-gas transition. The other end point 
is located at zero temperature (T = 0, n = 0.5), which 
is a common feature of the first-order PS described by 
DMFT'^'"^'*. The second-order charge ordering transi- 
tion line extends to zero temperature. Between this line 
(thin line in Fig. 4) and the dashed line lies an unstable 
area where the compressibility is negative. Metastable 
phases, either charge-disordered or charge ordered, ex- 
ist in the two patches between the unstable area and 
the boundaries. They arc intrinsic fcatmcs of first-order 
phase transitions and important for the properties of ma- 
terials near PS. In doped manganites, metastable phases 
have been observed through the hysteresis of magnetiza- 
tion and resistivity with respect to temperature"^^'^^, as 
well as through the resistivity relaxation phenomenon"^^. 
The study of metastable phases has provided valuable 
information about the PS in doped manganites. 
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FIG. 4. Phase diagram in the T — n piano for U = 2.0 and 
V = 1.4. The thin solid line is the second-order transition line. 
CO and CD denote charge-ordered and charge-disordered 
phase, respectively. PS denotes phase separation between 
charge-ordered and charge-disordered phases. The thick solid 
lines are the boundaries of the coexisting regime. In this 
regime, between the dashed line and the thin line lies a 
charge-ordered phase with negative compressibility. 



Electronic PS has been studied extensively in various 



models such as the Hubbard model, the t — J model, and 
the double-exchange model, etc. Most of the PS scenar- 
ios discussed so far rely on the magnetic exchange mech- 
anisms. Up to now, PS induced by pure Coulomb repul- 
sion is observed only in the spinless fermion modeP°. Our 
results show another example of PS caused by Coulomb 
repulsion only. In general, the Coulomb repulsion works 
against the PS, since a phase-separated state has a higher 
potential energy than the charge uniform state. Here we 
find that PS can also be driven by the charge-ordering 
transition which is in turn induced by pure Coulomb re- 
pulsion. When the charge orders to avoid the strong 
nearest-neighbor Coulomb repulsion in the high density 
regime, the Coulomb potential is strongly reduced which 
then allows for PS. It should also be noted that the PS 
described here exists mainly at finite temperatures, and 
comes from the effect of thermal fluctuations. In this 
respect, it is different from the ground-state PS of the 
spinless fermion model. 

In the following we discuss the effect of long-range (be- 
yond the nearest neighbor) Coulomb repulsion on the 
stability of PS between charge-disordered and charge- 
ordered phases. It is generally believed that long-range 
Coulomb repulsion will suppress a complete PS in the 
system, leading to an inhomogenious distribution of the 
electron density. Here we confine our discussion to the 
two-sublattice case and do not consider any incommen- 
surate ordering. Taking into account the long-range 
Coulomb repulsion, the total inter-site part of the Hamil- 
tonian (including the nearest-neighbor contribution) can 
be formulated as 



(15) 



In infinite dimensions, Eq. (15) reduces to its Hartree 
form. After performing appropriate scalings for the pa- 
rameters Vi^j, the effective mean-field Hamiltonian has 
the form: 



Hint =YniA {Vi {ub) + V2 {ha)) 
+ J2 ^"^) + ^2 (ns)) . 



(16) 



The effect of long range repulsion reduces to only the 
nerest-neighbor type and the next-nerest-neighbor type, 
which are represented by Vi and V2, respectively. If we 
denote (ua) = n + S and (ns) = n — S, then Eq. (16) 
further simplifies to 



iJint = ^ (Vl + V^2) n rii + ^ (y2 - VI) ^ UiA 



+ S{Vi-V2)J2niB. 



(17) 



This means that the long-range Coulomb repulsion in- 
troduces two effects. One is the enhancement of the av- 
erage repulsion as shown in the first term. The other is 
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the frustration effect caused by the next-nerest-neighbor 
type repulsion. In Eq.(17) this effect is reflected by the 

terms proportional to V2 — Vi. The first term can effec- 
tively suppress the coexisting regime, but cannot destroy 
PS at sufficiently low temperature^^. In the low tempera- 
ture Umit, the density difference of two coexisting phases 
reduces to zero, and the average repulsion will lose its 
effect on the PS. In contrast, the frustration effect in- 
duced by long-range Coulomb repulsion may destroy the 
PS completely, since PS crucially depends on the charge- 
ordering transition. In particular, if the frustration is so 
large that the effective nearest-neighbor V is less than Vc, 
as shown in the next section, there is no PS at all. There- 
fore we conclude that the long-range Coulomb repulsion 
may destroy the PS through its frustration effect. 

In the regime V > Vc, the charge order becomes robust 
with increasing V, as shown in Fig. I. Therefore, we 
expect that as V increases, the charge-ordered area in 
the T^n phase diagram will expand in both temperature 
and filling regime. However, the main shape of the phase 
boundary will remain unchanged. In particular, due to 
the close relation between PS and the charge ordering 
transition, the PS regime may well expand while keeping 
its main structure. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 



FIG. 5. The charge density difference between sublattice 
A and B as functions of T at (7 = 2.0, V — 1.4 for several 
electron fillings: n = 0.4 (dash), 0.45 (dot), 0.5 (solid), 0.65 
(dash-dot), and 0.7 (dash-dot-dot). 



Besides PS, we are also interested in the properties of 
the charge order in this regime. In Fig. 5, the order 
parameter \nA — is shown as a function of tempera- 
ture for different fillings. Except for n = 0.5, \nA — «b| 
shows a nonmonotonous behavior for all the fillings we 
studied. It first increases and then decreases upon lower- 
ing the temperature. In the filling regime 0.5 < n < 0.69 
where the ground state is charge ordered, \nA — ns] re- 
duces to a finite value at T = 0. For fiUings outside 
but close to this regime, the ground state is charge- 
disordered and a reentrant transition occurs at finite tem- 
perature. Compared with the reentrant transition, the 
non-monotonous change of the order parameter is a more 
general phonomenon. Experimentally, in systems with an 



ordered ground state but near the reentrant transition, 
such a non-monotonous change of the order parameter 
may well be observed. 



B. Reentrant chcirge ordering: V = 1.2 < Vc 

In this section, we discuss the case of V = 1.2, which 
is smaller than Vc ~ 1.32. For this repulsion strength, 
the ground state of the quarter-filled system is charge- 
disordered. As temperature decreases, the system shows 
a disorder-order-disorder type reentrant transition^ ^, as 
shown in Fig.l. When the electron filling moves away 
from n = 0.5, we find that the ground state is still 
charge-disordered, and that a reentrant charge ordering 
transition exists in an extended density regime 0.42 < 
n < 0.68. In Fig. 6, the curves for — versus T 
are shown for several fillings in this regime. It is seen 
that \nA —ns] varies non- monotonously as temperature 
decreases, and drops to zero at some finite temperature. 
Though the high temperature transition is continuous for 
all fillings, we find that the reentrant transition at lower 
temperatures has a different behavior for small and large 
fillings. For the filling n = 0.4 and n = 0.45, — Ub] 
changes smoothly at the reentrant transition tempera- 
ture, indicating that this transition is of second order. 
In contrast, in the curves for n = 0.55, 0.6 and 0.65, 
\nA — nsl drops abruptly near the temperature T = 0.05. 
Accompanied with this rapid change of the order param- 
eter an obvious slowing down of convergence is observed, 
which reduces the numerical precision significantly. A 
similar situation is found also for the filling n = 0.5. 
Since we have not found a hysteresis typical for first- 
order transitions, here we would not conclude that it is 
a first-order phase transition. More detailed studies are 
needed to elucidate this issue. 




FIG. 6. The charge density difference between sublattice 
A and B as functions of T at ?7 = 2.0, V = 1.2 for several 
electron fillings: n = 0.4 (dash), 0.45 (dot), 0.55 (solid), 0.6 
(dash-dot), and 0.65 (dash-dot-dot). 
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For both V = 1.4 and V = 1.2, the non-monotonous 
behavior of \nA — ns] versus T originates from the large 
spin entropy of the paramagnetic charge-ordered phase. 
For the charge-ordered phase, the two-fold spin degen- 
eracy on the occupied sites contributes a total entropy 
{N/2) In 2, while for the charge-disordered state near n = 
0.5, the entropy increases from zero with increasing T. So 
at low temperature, when the charge-ordered phase has 
a higher entropy than the disordered phase, the charge 
order develops as temperature increases. At higher tem- 
perature, where the entropy of the disordered phase ex- 
ceeds that of the charge-ordered phase, the charge order 
is reduced with increasing temperature. This is the main 
reason that leads to the non-monotonous behavior of the 
order parameter. Therefore, if the spin degeneracy of 
the charge-ordered phase is destroyed by the forth-order 
superexchange mechanism, it is doubted whether such a 
reentrant behavior is still present. Hellberg et al.^^ car- 
ried out a finite-temperature Lanczos study on a 4 x 4 
lattice. They found, however, that the reentrant tran- 
sition is stable when the superexchange effect is partly 
taken into account. 

In Fig. 7, the T — n phase diagram for ?7 = 2.0 and 

V = 1.2 is shown. The reentrant charge-ordering transi- 
tion exists in a regime from n w 0.42 to n « 0.68. Similar 
to the case oiV = 1.4, the charge ordering regime is also 
asymmetric around n = 0.5. Charge order appears only 
at finite temperatures. Here, one important difference 
from the V — lA diagram is that there is no PS near the 
reentrant transition line. For the fillings near the two ver- 
tical boundaries of the charge-ordered area, the order pa- 
rameter changes continuously to zero at the lower transi- 
tion temperature. This continuous transition is indicated 
by a thick line Fig. 7. In the regime 0.5 < n < 0.65, 
the high temperature transition is continuous, while at 
low temperatures, the order parameter disappears more 
abruptly. In the phase diagram we schematically denote 
such an abrupt transition by a dotted horizontal line. 
Due to the severe critical slowing down of convergence 
at this transition, it is difficult to obtain the transition 
temperature accurately. 

Comparing Figs. 4 and 7, we see that the T — n phase 
diagrams for V = 1.4 and V = 1.2 are topologically dif- 
ferent. Near the critical Vc, the T — n phase diagram, 
including the stability of PS and the reentrant transi- 
tion, is very sensitive to V. When V increases from 

V < Vc, both the n and T regime of the charge-ordered 
area expands, and the reentrant transition temperature 
decreases. But the ground state remains disordered, as 
shown in Fig. 7. At V = Vc, the ground state at n = 0.5 
is expected to first turn charge-ordered. When V is even 
larger, the filling regime of the charge-ordered ground 
state extends towards larger n, while keeping its left end- 
point n = 0.5 unchanged. At the same time, a PS area 
where charge-disordered and charge-ordered phases co- 
exist emerges near the reentrant transition line in the 
n < 0.5 regime. In this way, the phase diagram shown in 
Fig. 7 for small V evolves into that in Fig. 4 for large V. 



Phenomena such as charge ordering, reentrant transi- 
tion, and PS have been observed experimentally in doped 
manganites. The mechanisms for these phenomena are 
more complicated and are topics of intensive research. 
The interesting point here is that starting from an elec- 
tronic model that only takes into account the on-site and 
nearest-neighbor Coulomb repulsion, we are able to ob- 
tain rich phase diagrams that include all these phenom- 
ena. Previous studies^^ indicated that the occurence of 
a reentrant charge ordering transition does not depend 
on details of the localization mechanism. Here, when 
compared with other first-order phase transitions stud- 
ied before, we find that within DMFT, the structure of 
first-order transitions and their temperature behavior are 
very similar. They are also independent of the specific 
mechanism. It would be interesting to do such a com- 
parison for those real materials where different first-order 
phase transitions occur. 
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FIG. 7. Phase diagram in T — n plane for U = 2.0 and 
V = 1.2. The solid line is the second-order transition line. CO 
and CD denote charge-ordered and charge-disordered phase, 
respectively. The dotted line denotes an abrupt change of the 
order parameter. 



IV. SUMMARY 

In this paper, we studied the paramagnetic phase of the 
single band extended Hubbard model near quarter filling. 
In the framework of DMFT, the effective Anderson im- 
purity model is solved using the exact diagonalization 
technique. Based on the previous T — V phase diagram 
for n = 1/2 (Ref."'^^), we concentrate on the phase dia- 
gram in T— n planes for two specific value of the inter-site 
repulsion V: V = 1.4 and V = 1.2. located on the two 
sides of the zero temperature critical value ~ 1-32. In 
both cases, charge order exists in an extended regime of 
filling near n — 1/2. Except at some special points (such 
as V = 1.4, n = 1/2), the order parameter of charge or- 
der \nA — ns] changes non- monotonously with tempera- 
ture. For V = 1.4, in the filling regime 0.39 < n < 0.5 
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and near the reentrant transition temperature, we find 
PS between a charge-disordered and a charge-ordered 
phase. Information on this PS, including the existence 
of a metastable state and the first-order transition hue, 
is described by the "S" -shaped structure in the contin- 
uous n — /i curve. At zero temperature, with increas- 
ing n, the charge-disordered ground state changes into 
a charge-ordered one abruptly at n = 1/2. Our analy- 
sis suggests that long-range Coulomb repulsion may de- 
stroy this PS through its frustration effect. For V = 1.2, 
the ground state is charge- disordered for all fillings. The 
reentrant charge ordering transition is observed in the 
regime 0.42 < n < 0.68. It becomes rather abrupt in the 
regime 0.5 < n < 0.65. Relevance of our results to the 
doped manganites is discussed. 
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